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(<-^ , Abstract. Understanding of the so-called 2 particle-2 hole (2p-2li) effect is an urgent program in neutrino interaction physics 
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for current and future oscillation experiments. Such processes are believed to be responsible for the event excesses observed 
by recent neutrino experiments. The 2p-2h effect is dominated by the meson exchange current (MEC), and is accompanied 
by a 2-nucleon emission from the primary vertex, instead of a single nucleon emission from the charged-current quasi-elastic 
(CCQE) interaction. 

Current and future high resolution experiments can potentially nail down this effect. For this reason, there are world wide 
efforts to model and implement this process in neutrino interaction simulations. In these proceedings, I would like to describe 
how this channel is modeled in neutrino interaction generators. 
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O : MESON EXCHANGE CURRENT (MEC) 

C I The MEC is conceived to be responsible for the so-called "dip region". In electron scattering (for review, see [1]), 

inclusive electron scattering measures the out-going electron from the interaction. Since the beam has well-defined 
ff^ I energy, an energy transfer co and a 3-momentum transfer q are reconstructed with high accuracy. The inclusive cross 

K^ ■ section is the differential cross section as a function of energy transfer In this space, with moderate electron beam 

\l ' energy, one can see 2 bumps, representing quasi-elastic (QE) scattering peak and A-resonance peak. These are often 

reasonably modeled, however, many theories fail to reproduce a dip between these 2 bumps. This is the dip region. 
By adding the MEC, some models successfully reproduce inclusive cross section data including dip region [2]. The 
contribution of MEC is larger in the transverse response than the longitudinal response. 

Therefore, the importance of MEC is known from electron scattering data. MEC is an interaction involved in 2 
nucleons, or 2-body cuiTent, and it is classified in "2 particle-2 hole (2p-2h)" effect. Here, a weak boson from the 
leptonic current is exchanged by a pair of nucleons (2-body current), and believed to lead to 2-nucleon emission. 
The importance of this process in neutrino interactions was first pointed out by Martini et al. [3] shortly after 
the MiniBooNE experiment showed their CCQE double differential cross section [4]. Several groups successfully 
reproduces the MiniBooNE CCQE or neutral current elastic (NCEL) cross section data [5] by adding the MEC in their 

^ . models [6, 7, 8, 9, 10]. However other approaches appear to work equally well [11, 12]. 

^.' To understand weather the MEC is the responsible process for enhancement of recent neutrino data [4, 5, 13, 14, 15], 

we need to compare more kinematic details of the data with simulations. Here, we desciibe how this process is 
implemented in GENIE neutrino simulation generator [16]. 

Through this note, I assume charged-cuiTent interaction on carbon target, also, a neutiino scatters with a neutron- 
proton pair (hence outgoing particles from the primary vertex are one muon and two protons). 

MODELING MEC NEUTRINO INTERACTION IN GENIE 

Figure 1 shows a cartoon of how to model MEC in GENIE. Any interaction can be modeled with these three parts. 

• Leptonic model 

• Hadronic model 

• Final state interaction (FSI) model 



1 . Leptonic model 
(Dytman model) 

V 




2. Hadronic model 
{Nucleon cluster model) 



3. FSI model 
(hA model) 



FIGURE 1. Basic strategy of modeling MEC in GENIE. 



The leptonic model defines the differential cross section of the outgoing lepton for a given energy-momentum 
transfer So this part completely specifies the lepton kinematics. There are several models available, and a neutrino 
interaction generator may choose one. Then using this energy-momentum transfer, the physics of the hadronic system 
needs to be specified. The problem is that there is no microscopic hadronic model of MEC is available. Distribution of 
initial nucleon kinematics, presence of correlations, type of nucleon pairs, etc, are key ingredients producing outgoing 
nucleons. This needs to be modeled. Finally, all outgoing particles experience FSI in the target nuclei. Observables 
always depend upon FSI, unfortunately. 



Leptonic model, Dytman model 



To specify leptonic kinematics, GENIE collaboration developed the Dytman model. The Dytman model is motivated 
by the Lightbody model [17]. In the Lightbody model, QE and A peaks are given by Gaussian distributions, then the 
MEC part is modeled as a Cauchy distribution between those 2 distributions. In the Dytman model, MEC is modeled 
as a Gaussian distribution between 2 peaks, QE and A, to fill the dip region of electron scattering inclusive cros s section 
data. The target nuclei dependence is linearly given. 

The differential cross sections for electron and neutrino scattering of the Dytman model are made using only the 
Sachs magnetic form factor (G^ p/„), to emphasize the transverse nature of MEC [8]. The Sachs magnetic form factor 
has power 6 (dipole form factor is power 2), motivated by the deuteron elastic scattering form factor Here, electron- 

nucleon scattering can be written (where T = '^i)^ 
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and the neutrino-nucleon CC cross section can be written. 
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In tMs way, one can universally treat MEC in the electron scattering and the neutrino scattering processes. 

The goal is to tune the MEC strength from the electron scattering data, and apply the same model to the neutrino 
cross section prediction. The agreement is tested with electron scattering data. Figure 2 (left) shows an example of 
inclusive cross section prediction. The data points are taken from Ref. [18]. The model tuning using electron scattering 
data is ongoing. In this moment (GENIE v2.7. 1), the strength is chosen so that it agrees with MiniBooNE and NOMAD 
data. Since NOMAD [19] does not observe the enhancement, the Dytman model is chosen to turn off monotonically 
from 1 GeV to 5 GeV (Fig. 2, right). 
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FIGURE 2. Dytman model predictions. The left plot shows data-simulation comparison of the 560 MeV electron inclusive cross 
section with 60 degree scattering. The MEC contribution can be seen between QE and resonance peaks. The right plot shows energy 
dependence of the neutrino total cross section with a function of neutrino energy. 






FIGURE 3. Nucleon cluster model in GENIE. First, 2 nucleons are chosen from the Fermi sea (left). Then 2 nucleons and energy- 
momentum transfer 4-vector make the hadronic system (middle). Finally, the system is boosted back, and 2 outgoing nucleons are 
generated in the lab frame (right). 



Hadronic model, nucleon cluster model 



For the hadronic model, the GENIE collaboration developed the nucleon cluster model. The hadronic model is 
important because this specifies outgoing nucleon types and kinematics. The nucleon cluster model is a naive approach, 
but reasonable (a similar model was independently developed in NuWro [20]). The procedure is depicted in Figure 3. 
Here 2 nucleons are independently selected from the Fermi sea, and the sum of them makes a "nucleon cluster" (Fig. 3, 
left) with momentum P where the leftover nuclear recoils with momentum —P. Here, no Pauli blocking is applied. 
Also no separation energy is considered. Then, this nucleon cluster and energy-momentum transfer 4-vector form the 
center-of-mass system, called the "hadronic system" (Fig. 3, middle), and isotropic decay in the hadronic system is 
boosted back to the lab frame to produce 2 outgoing nucleons (Fig. 3, right). 
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FIGURE 4. Comparisons of GENIE and MiniBooNE CCQE double differential cross section. Here all plots are a function of 
muon kinetic energy, in a slice of muon scattering angle. CCQE is defined "one muon, no pion, any number of nucleons" in the 
detector. From left to right, 0.8 < cosdfi < 0.9, 0.5 < cosOfj, < 0.6, -0. 1 < cosd^ < 0.0. 

Final state interaction (FSI), hA model 

For the FSI model, the MEC model in GENIE uses the hA model as a default. More detail can be found else- 
where [21]. 



Results 

Here, simulation results are generated by folding with the MiniBooNE v^ flux file [4]. 

Lepton kinematics 

Figure 4 shows the comparison of GENIE with the MiniBooNE CCQE double differential cross section data. In 
MiniBooNE, cross section channels are defined from the final state particles in the detector instead of interaction types 
at the primary neutrino interaction vertices [22]. Here MiniBooNE defines CCQE as "one muon, no pion, any number 
of nucleons" in final state particles in the detector. The Dytman model has a good agreement with the MiniBooNE data 
at a moderate scattering angle, but it slightly overestimates at forward going muons, and underestimates back scattered 
muons. The tuning with larger set of electron scattering data is ongoing. 



Hadron kinematics 

It may be more interesting to take a look hadronic system prediction of GENIE. Figure 5 shows examples of the 
GENIE prediction of nucleon kinematic distributions. Here histograms are arbitrarily normalized to emphasize the 
shape differences. The left plot shows the leading nucleon momentum. Since MEC shares energy-momentum transfer 
with 2 nucleons, the maximum available nucleon momentum is lower than what simple CCQE predicts. This indicates 
that it may be a little harder to see the tracks of nucleons from the MEC. The middle plot shows the opening angle of the 
2 leading nucleons. Although it is isotropic for the CCQE and resonance processes where the second leading nucleon 
is created by FSIs, MEC clearly predicts larger angles at the MiniBooNE neutrino beam energy (^ 800 MeV) due to 
the nucleon cluster model. The right plot shows the sum of all nucleon kinetic energies. Such quantity may be more 
useful if experiments cannot identify multiple nucleon tracks [20]. Here, GENIE MEC model makes a characteristic 
peak around ^Q3 GeV due to the Dytman model. This may be compared with predictions from other generators based 
on different leptonic models [20] . 
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FIGURE 5. Predictions of GENIE MEC model on nucleon kinematics. The left plot shows the momentum of the leading nucleon, 
the middle plot shows the opening angle of the 2 leading nucleons, and the right plot shows the sum of all nucleon kinetic energies. 
Note histograms are arbitrarily normalized to emphasize the shape differences. 

OTHER NEUTRINO INTERACTION GENERATORS 

The model of the 2p-2h effect is also developed in NuWro [25] and GiBUU [26] neutrino interaction generators. 



NuWro 

NuWro predictions of the nucleon kinematic distributions are shown in Ref. [20]. Hadronic kinematics predictions 
are related with not only the hadronic model, but also the leptonic model, because the leptonic model specifies the 
energy-momentum transfer given to the hadronic system to emit nucleons. Currently, NuWro prepares 3 different 
leptonic models, np-nh model based on Marteau et al. [27], Valencia model from Nieves et al. [7], and Bodek et aZ.'s 
transverse enhancement model (TEM) [8]. 



GiBUU 

GiBUU [9] has more emphasis on the FSI part, by taking account the phase space evolution. Interestingly, GiBUU 
predicts no MEC contribution in "one muon, zero pion, zero proton" CC sample. GENIE agrees with this result. This 
opens a new way to measure genuine CCQE interactions, since precise counting of outgoing nucleons may nail down 
the primary process. NOMAD [19], K2K [23], and SciBooNE [24] showed the zero and one proton counting from 
their CC samples, but future high resolution experiments, such as MicroBooNE [28], can push this idea further. 



CONCLUSION 



Table 1 shows a comparison of MEC models in various generators. Clearly, the hadronic model is the key to find the 
MEC processes in the neutrino scattering program. Are there any correlations in initial nucleon kinematics? What 
kind of pairs are favored, proton-neutron pairs or neutron- neutron pairs? and by how much? How to share the energy- 
momentum transfer between 2 nucleons? etc. All these are open questions! Theorists are encouraged to provide this 
information to experimentalists. 

The GENIE collaboration is preparing to release the next frozen version, v.2.8.0. Here, the MEC model described 
in this paper will be included. There is an ongoing effort for the further tuning of this channel based on the larger set 
of electron scattering data. 



TABLE 1. Comparison of MEC models in neutrino interaction generators. 



GENIE 



NuWro 



GIBUU 



Leptonic model 

Hadronic model 

initial nucleon momentum 

initial nucleon momentum correlation 

initial nucleon spatial correlation 

initial nucleon pair 

FSI model 



Dytman model 

nucleon cluster 

Fermi sea 

none 

none 

n-p:n-n=l:4 

isospin ansatz 

hA model 



TEM, np-nh model, 

and Valencia model 

nucleon cluster 

Fermi sea 

none 

none 

n-p:n-n=9:l 

short range correlation 

cascade model 



Transverse projector 

phase space density 

Fermi sea 

none 

2 nucleons are generated 

at the same location 

n-p:n-n=12:5 

statistical average 

BUU transport 
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